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Abstract 

We examine in detail the cosmology based on quantal (Bohmian) trajectories as suggested in a recent 
study [1]. We disagree with the conclusions regarding predicting the value of the cosmological constant 
A and evading the big bang singularity. Furthermore, we show that the approach of using a quantum 
corrected Raychaudhuri equation (QRE), as suggested in [1], is unsatisfactory, because, essentially, it uses 
the Raychaudhuri equation, which is a kinematical equation, in order to predict dynamics. In addition, 
even within this inconsistent framework, the authors have adopted unjustified assumptions and carried 
out incorrect steps leading to doubtful conclusions. 


It has been widely known that the cosmos, at a large scale, is to be homogenous, isotropic and spatially 
flat and thus could be described by Friedmann-Robertson-Walker (FRW) metric whose line element in 
the comoving cartesian coordinates is given by, 

dx’' = —+ a(t)^ [(dx^)^ + , (1) 

where x^ is the four dimensional coordinate, x^ = (a;° = ct, x^ = x, x'^ = y, x^ = z), and a{t) is the 
scale factor. 

The metric tensor can be easily read from Eq.(l) to be diagonal and given by, 

gfj,^ = Diag [-1, a{tf, a{tf, a{tf] . (2) 

The scale factor a{t) can be determined by applying field equations of General Relativity (GR), namely: 
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where and R are the Ricci tensor and scalar respectively. As to the energy-momentum tensor 
which in our case represents a perfect fluid having density p, pressure p and velocity t/^, it assumes the 
form 

{p+ ^^UyUv+pgf,^. ( 4 ) 

It is worth mentioning that we are using a metric with signature (—, -I-, -I-, -I-), in order to be the same one 
used in [1]. We also do not set c, the speed of light, equal to unity for the sake of clarity and to keep all 
fully dimensionfull factors apparent, like the Newton’s gravitational constant G. For later manipulations, 
we recall that the Hubble parameter H is defined to be d/a. 

The current observations [2]-[5] ranging from Type lA supernova observations, cosmic microwave 
background radiation (GMBR ) data and baryon acoustic oscillations indicate that the matter-energy 
content of our universe is comprised of dark energy filling 72% , probably in the form of a cosmological 
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constant A whose physical origin is debatable, dark matter occupying about 23% whose physical non- 
baryonic nature is so far not determined, and the rest as an ordinary baryonic matter which could be 
luminous or not. 

The authors in [1] used a quantum corrected Raychaudhuri equation (QRE) to derive a modihed 
Friedmann equation which enabled them to extract a value for the cosmological constant A as well as to 
evade the big-bang singularity. The QRE was obtained by replacing geodesics with quantal (Bohmian) 
trajectories. In order to show why, to us, this procedure is invalid, it is appropriate to present briefly the 
essential ingredients for Raychaudhuri equation (RE) as can be found in many books on general relativity 
such as [6, 7]. 

The RE is concerned as regards how a congruence of curves might evolve. The evolution can be 
attributed to a change of the quantity which is defined as, 

B^j,v = (5) 


where is a tangent vector along the curve parameterized by affine parameter r along the curve. For 
convenience, we suppress the dependence on the other parameter expressing how far the neighboring 
curves are apart from each other. The rate of change of along the curve is given as, 

= V^DxB^, = D, {V^DxV^) - {D,V^) (DxV^) - , (6) 


where Di^ denotes the covariant derivatives and RafiXv is the Riemann curvature tensor * Having chosen 
a congruence of timelike geodesics we get 


Defining 6 = B^^, 


DB/j^i 

Dt 


= -B\B^x-Ra^.XuV^Vr 


and plugging it into Eq.(7), we get 
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Now we decompose B^^, into its irreducible parts as, 


B^li — 


lep, 


fxi/ ~r 1 


(9) 


where the projector P^i, projecting onto the 3-dimensional subspace orthogonal to is defined as 


P^^v — gfiv + VfjYu. ( 10 ) 

The three terms in Eq. (9) have distinct physical interpretations in that the trace part 9 describes 
expansion, the symmetric traceless part = i {B^y + Bi^^) — \9P^^ describes shear and the skew 
symmetric part ^ — B^^) describes rotation. 

When plugging Eq.(9) into Eq.(8) we get the celebrated RE: 


D/9 1 

- 02 _ ^ ^ _ R ( 11 ) 
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RE as it stands is a kinematic equation devoid of any dynamics. In order to get in touch with the 
physical world, the tensor R^i, should be supplied as the one that satisfies GR equations and the specihed 
congruence must be that of timelike or null geodesics. 

In a standard cosmology, we have FRW metric in the comoving frame where the cosmic fluid can be 
considered to be at rest i.e. = (c, 0,0, 0). The quantities and uj^i, are vanishing as can be easily 
verified, while the expansion parameter 6 turns out to be 3i7. Upon using GR equations (Eq.(3)) for i?oo 
we get, 

7?oo = ^4 + 3p) , (12) 

*The convention followed for the Riemann curvature and Ricci tensors are respectively = {ppXZp ~ 

(di^rjp -|-rj„r“p) and 
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Plugging the result for Rqq in Eq.(ll) we get the ‘Second’ Friedmann equation: 

47rG 


H = -H^ - 




( 13 ) 
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Expressing the FRW Ricci scalar as R = [ — I —^ ) ^.nd writing i?oo from (Eqs. 3 and 4) as 


pc^, we get the ‘First’ Friedmann equation: 


= 
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(14) 


To summarize, we can get the two basic Friedmann equations of cosmology, Eqs.(13,14), through the 
RE (Eq.(ll)). However, this occurs only after imposing the GR field equations and using a congruence of 
timelike geodesics that coincide with the motion of the cosmic fluid. In this specific case, the expansion 
parameter 9 is obliged to be 3H. Thus, any obtained equation governing H should be consistent with the 
GR field equations. Having chosen different congruences, not necessarily geodesics, then the expansion 
parameter 9 would not be 3H. Hence, the derived equation for 9 based on Eq.( 6 ) should be solved to get 
the expression of 9. It is important to stress that one can not derive new dynamics from RE other than 
that contained in the GR field equations. 

The authors of [1] have used the idea of quantal Bohmian trajectory in conjunction with RE. Although 
the mere idea of quantal Bohmian trajectory is a problematic one, as we shall show later, but nevertheless 
let us continue with our understanding of the authors’ arguments. The authors started with a massive 
scalar field $ coupled non-minimally to gravity. The equation of motion for $ is taken to be in the form. 


□ - 


F 2 


eR $ = 0, n = g^''D^D,, 


(15) 


To get the velocity field associated with a Bohmian trajectory, one can use the polar decomposition of $ 
as. 


then Eq.(15) would imply 


$ = Re R 


nv 

gf^’'iD^S){D,S) = -m^c^ + re— + hhR, 


2gf^’' {D^TZ){D^S) = -TZDS. 

The velocity field associated with the quantal (Bohmian) trajectory is defined to be: 


and as a consequence of the first of Eqs.(17) we obtain, 


U^U„, = -c" 


DR 


m 


—y eR. 


(16) 
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(19) 


It is clear that the velocity field U^, as its norm is defined in Eq.(19), is not guaranteed to be always of a 

DR 

timelike type. It could flip into a spacelike vector due to the presence of the terms —tt —— and — 5 - e R. 

m‘‘ R ^ 

This is a serious drawback that might turn the whole approach into being unphysical. Putting aside this 
difficulty, let us continue with the usual procedure and compute the quantity [DfJJi,) U'^ to find: 


{D^U,) t/" = 


2 m 2 ^ 


f un 




( 20 ) 


Using Eq.(20) and Eq.( 6 ), one can obtain the QRE (the modified form of the RE applicable to these 
non-geodesic quantal trajectories) as: 


D9 

Dt 




m 


+ —^eDR + A{9,h,R,n, U>^), 


2 m 2 


A{9,h,R,n,U>^) = - 


3 UM’' 2m? 


eU^D^R -k 


f un 


( 21 ) 
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dx^ 

where r is the parameter defined by = ——. Let us compare the QRE of (Eq. 21) with the 

UT 

corresponding QRE in [1]: 

Dt 6 m‘‘ \ TZ J 

We find that the two QRE equations do not coincide in that the and the A terms in (Eq. 21) are 
absent in (Eq. 22), and the numerical factors appearing in front of are also different even though the 
precise value of these factors are not relevant for our intended considerations. We trace the absence of 
the A—term in (Eq. 22) to the ‘improper’ definition in [1] of the trace part 

9 = (23) 


in that it uses the projector as defined in (eq. 10) instead of using the FRW- metric for the 
correct definition of the trace part: 

9 = (24) 

This explains also the form of the two h^-terms in (Eq. 22) and why they differ from the corresponding 
terms in (Eq. 21). We note here, considering the new velocity squared-norm (Eq. 19), that the definition 
of the projector should be amended to be: 

Pfiv = Qfiv JJaJJ ( 2 ^) 


We note however that using either or g^^" would amount in [1] to the same result for the h^-terms 
in Eqs (21, 22) upon neglecting higher orders of h. This comes because the velocity in [1] was defined as 


= h 


P^(|) 


= h 


Df^S 


and so the U^U^-term in the ‘imprecise’ definition of the projector (Eq. 10), 


m m 

in [ 1 ], was considered as a quantum correction producing a higher order term in h (this argument would 
not apply for the A-term since the P^^’'- contribution vanishes and one needs to compute explicitly the 
C/'^P'^-contribution in order to evaluate the full g'^^-contribution). 

As to the term containing describing shear, one should not discard it for the specific choice of 
the velocity field with norm given by Eq.(19). Actually the non-zero components for turn out to 
be dependent on both H and S. When assuming S depending only on t in order to be consistent with 
homogeneity and isotropy, then the non-zero components for can be easily computed and we get: 


o'oo 


S 





CTii = CT22 = 0-33 


So? 

3mc^ 


(26) 


Again, putting aside the discrepancies in the two QREs, the authors of [1] have later manipulated 
their QRE (Eq. 22) confusingly, to us, several times: 


First: They have identified 9 with SH which is inexact, in our opinion, because the particle motion 
described by the velocity field U^, as given in Eq.(19), is not a geodesic motion. Assuming S to be t- 

3 • S 

dependent in accordance with homogeneity and isotropy, then 9 is given hy 9 = - - H S - 

Thus, the identification of 9 with 3iL and interpreting the resulting equation for H as a, new 
Friedmann equation is not accurate. 


Second: They have replaced RaxU^V^ by (pc? 3p) which is not completely true as it should be 

c 

AttG 

replaced by — 4 —(pc^ + 3p)—^ for a t-dependent S. 


Third: They have assumed a value for TZ equal to exp where r = yxc^ + y^-l-^, and they 

considered Lq as the present size of the universe. The position dependence of TZ is not in accordance 
with homogeneity and isotropy of the FRW-metric, as the latter two features admit only time 
dependence for any scalar function. 
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Fourth: They have identified the contribution 






mz\ 


) as a cosmological constant whose value 


would be Aq = 


out to be, 


2m?c'^ y TZ J 


2m? \ n J 

which is not a constant as stated in [1]. The value of Aq turns 




2m‘^a^Ly 


6 + ^ (2r^ - 3Ll) 


(27) 


Fifth: They have carried out the calculations with ordinary derivatives instead of covariant derivatives, 
especially in computing Aq which explains how they got a diffrent result {^)- Apparently, the 
authors in [1, 8] considered using the spacetime metric as a fixed background which would allow 
them to use ordinary derivatives instead of covariant ones. We think this is not correct because one 
will lose the general covariance. 


Leaving aside all the above mentioned debatable points, which are enough to put into doubt the work 
carried out in [1], let us, for the sake of complete judgement, continue the discussion of the QRE obtained 
in [1]. After having identified 6 with 3iL, RaxU^V^ with + 3p), and after having dropped the 

and A-terms in (Eq. 21), and after having discarded the term which was identified as a cosmological 
constant (-fortunately not relevant for the foregoing discussion-) we get, assuming the validity of the GR 
field equations and using the equation of state p = w pc?, the following equation: 

H = -\[l+u^)H^ + ^cUR. (28) 

Apart from an irrelevant, for our discussion, numerical factor in front of f?, Eq.(28) coincides with the 
corresponding equation in [1]. Now, regardless of any quantum theory of gravitation having GR as a 
zeroth-order approximation, one can substitute, up to a multiplicative numerical factor, T -the trace of 
the energy-momentum tensor- for R in Eq. (28) as the difference would constitute a higher order in h. 
Therefore, one can guess certain basic features that should be present in any contribution proportional 
to DR, namely: 

First: It should vanish for the case of radiation where (w = |) since T is zero for radiation. 

Second: It should vanish for the case of dark energy (cosmological constant where w = —1) since T is 
constant in this case. 


Third: It should vanish for the case of w = — ^ that corresponds to the case of Milne universe which is 
equivalent to an empty universe dominated by negative curvature term. 

SttG- 

However, assuming the GR field equations (which lead to R = - -r—T where T = = 

—(I — 3a;)), and the continuity equation (p = —3iL (I-I-w) p) and the Friedmann equation (iL = 
— I (I -I- w) H^), one can evaluate the contribution proportional to DR in Eq. (28) to get: 

H = ^ e (l - 90;^ (I + a;) H\ (29) 

which disagrees with what was obtained in [I]: 


H — -- (I +uj)H'^ -^ e (I + w) 

2 m^m 


6(1 


9 81 

- w) —— (I -I- w) -I- 18 




(30) 


The result in Eq.(29) possesses all the required features while that of Eq.(30) does only posses the 
second feature. The correct equation (Eq.29) reduces exactly to the GR equation in case of radiation, 
which is the most relevant case for the early universe. Having no deviation from GR in case of radiation- 
dominated universe implies that if we trace the cosmic evolution backward in time we would hit the 
big-bang singularity in a finite time. This finding does not agree with the conclusion derived in [I] based 
on the equation (30) where the big-bang singularity is completely evaded. 
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Avoiding the big-bang singularity can be simply explained through interpreting the equations govern¬ 
ing H, like Eq.(29) or Eq.(30), as a dynamical system. The fixed points associated with the dynamical 
system, where H = 0, can determine qualitatively the dynamical behaviour of H. For more elaborations 
on the notion of dynamical system and fixed points in a one dimensional flow, one can consult [9]. 

Creating attractive fixed points in the past at a finite value of H prevents H from running away to 
infinity in a finite time (a situation characterizing a big bang singularity). The system, instead, will 
be directed towards the hxed point in an infinite time lapse. The equation (29) ceases to create these 
attractive fixed point in the past for the case of radiation dominated universe (w = ^). The sign of the 
contribution proportional to in Eq.(29), responsible for creating non-trivial fixed points, is controlled 
by the factor e (l — 9 (1 -|- oj). This factor is zero for w = ^ or — 1, which are the most relevant cases 

for the early universe. 

In more details, one can examine, for positive epsilon, the relative sign of both contributions that are 
proportional to and H'^. The interesting values for tu should lie in the interval [—1, which contain 
the mostly motivated physical cases. When w s] — i, i[, the terms proportional to and H'^ have the 
same sign and thus no possible fixed points other than H = 0 can be found. On the other hand, the 
opposite is true for to €] — 1, — where the two contributions are against each other and a fixed point 
other than H = 0 can be generated. Regarding the case of negative e, the consequences concerning the 
two intervals a; g] — and ui g] — 1, — are swapped with each other. To sum it up, we did not find 
a deviation from GR when a; = i or — 1, and these two cases are very relevant for the early universe. In 
less relevant cases for the early universe, where uj ^ ^ or — 1, one can have a deviation from GR which 
is not so significant. 

A serious persistent drawback for all values of w, linked to the velocity field with norm defined in 
Eq.(19), is that it can not be guaranteed to be always timelike. The value of the squared-norm of in 
the case TZ = exp (—r^/Lp) is given by 


U^^U^ = -c2 -p 


rm? 


/ 4,^2 _ 02^2 \ 

V J 




(I -3w)ij2, 


(31) 


It is clear, now, that it is difficult to tune the vector type to be timelike through the whole history and 
in different patches of the universe. It suffices to take a large value for r in order to enforce to be 
spacelike. 

We conclude that the work carried out in [1] is questionable, particularly over the following two points: 


• Extracting dynamics from a kinematical relation like RE in conjunction with an unlikely identifi¬ 
cation of the parameter 9 and suspect calculations leading to problematic conclusions. 


• The velocity field associated with the quantal (Bohmian) trajectory can not be guaranteed to 
be always timelike, and this difficulty alone is sufficient to make the whole approach unphysical. 


Hence, in our opinion, the conclusions presented in [1], concerning evading the big-bang singularity and 
predicting a value for the cosmological constant A, are disputable. 
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